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fJJ 1 ABSTRACT 



Aims. We investigate conditions under which kink magnetohydrodynamic waves propagating along photospheric uniformly twisted 
flux tubes with axial mass flows become unstable as a consequence of the Kelvin-Helmholtz instability. 

Methods. We employed the dispersion relations of kink waves derived from the linearised magnetohydrodynamic equations. We 
assumed real wave numbers and complex angular wave frequencies, namely complex wave phase velocities. The dispersion relations 
were solved numerically at fixed input parameters and several mass flow velocities. 

Results. We show that the stability of the waves depends upon four parameters, the density contrast between the flux tube and its 
environment, the ratio of the background magnetic fields in the two media, the twist of the magnetic field lines inside the tube, and 
(— I _ the value of the Alfven-Mach number (the ratio of the jet velocity to Alfven speed inside the flux tube). We assume that the azimuthal 

component of the magnetic field in the tube is proportional to the distance from the tube axis and that the tube is only weakly twisted 
(i.e., the ratio of the azimuthal and axial components of the magnetic field is low). At certain densities and magnetic field twists, 
an instability of the Kelvin-Helmholtz type of kink (m = 1) mode can arise if the Alfven-Mach number exceeds a critical value. 
In particular, for an isolated twisted magnetic flux tube (magnetically free environment) at a density contrast (the ratio of the mass 
CZ) . density of the environment to that of the tube itself) equal to 2 and a magnetic field twist (defined as the ratio of azimuthal magnetic 

field component at the inner surface of the tube to the background magnetic field strength) equal to 0.4, the threshold Alfven-Mach 
number has a magnitude of 1.250075, which means that for an Alfven speed inside the tube of 10 km s _I the jet velocity should be 
higher than 12.5 km s to ensure the onset of the Kelvin-Helmholtz instability of the kink (in = 1) mode. Speeds of that order can be 
detected in photospheric tubes. 

Conclusions. The observed mass flows may trigger the Kelvin-Helmholtz instability of the kink (m = 1) mode in weakly twisted 
photospheric magnetic flux tubes at critical Alfven-Mach numbers lower that those in untwisted tubes if the magnetic field twist lies 
in the range 0.36-0.4 and the flow speed exceeds a critical value. A weak external magnetic field (with a ratio to the magnetic field 
inside the tube in the range 0. 1-0.5) slightly increases that critical value. 
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1 . Introduction l2008l Wa llace e t al. I20TU1 Madiarska l20TTl K amio et al. I2TJTT1 

Tian et al. EOTTl Srivastava & Murawski lSOTTT l. 

An additional feature of the flux tube is its twist. There 

With the wealth of recent high-resolution observations from is observational evidence for a twist in the solar atmosphere. 

b ■ SDO, Hinode, Stereo, TRACE, SOHO, and RHESSI as well as Aschwanden et al . 47999) confirmed the existence of twisted 

from ground-based observations, the concepts of the magnetic loQps in ±e corona by detecting a post . flare i oop system in ex- 

flux tube and coronal loops have become increasingly impor- . , . , s, . , „ , ~ , 

c , ,. f . , 1 , J, treme ultraviolet wavelengths (171 A) with the Coronal Explorer 

tant for understanding explosive phenomena such as solar flares ™„ . „„ x „ , ,. , , . , , , . 

, . • jfr.i 11 (TRACE). Rotational movement along a loop as observed in 

and eruptive prominences and also for the problem of the so- \ A c ^ i_ * i_ i_ m. 

, ■ • ■ • „ , • , • , the active region NOAA 8608 of the photosphere by Chae 

lar corona heating mechanism. The observations reveal a wide , i „n».i . ,. . 5 ■ • 

. , , . . , „ . et al. (2001) indicates that there is a twist in the magnetic 

range of shapes and sizes, ranging from the small-size mag- „ , , ,. r. „ . . , i~rrnxi i.- 

F „ . J-.-, , , field lines. More recently, Srivastava et al. (2010), using multi- 

netic flux tubes found in the granular network, sunspots, coro- , , .. J . . '' . 5 n_i 

. , b , , , wavelength observations of SOHO/MDI, SOT-Hinode/blue- 

nal loops, and arcades with many loops to huge loops such as . , 5 , „ , , = , _ .. „„ ' = 

prominences. The registered mass flows speeds in these mag- ^f™ Jfi* , A >' ban , d ( 430 u 5 A >' CaI ! ( 3968 ^ 

netic tubes range from a few tens of kilometers per second in and TRACE 171 A, presented the observational signature of a 

the photosphere/chromosphere (Shibata et al. [20071 Katsukawa T^7J"»ft~ loop in AR 10960 during the period 

et al. [20071 De Pontieu et al. [2003 up to hundreds or thou- 04:43 UT-04:52 UT on June 4, 2007. 

sands of kilometers per second in coronal holes and X-ray jets It is well-known that the hydromagnetic flows are generally 

(Shimojo & Shibata 120001 Cirtain et al. 120071 Ofman & Wang unstable against the Kelvin-Helmholtz instability when the flow 
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speed exceeds a certain threshold/critical value (Chandrasekhar 
119611 1. Recent observations of Kelvin-Helmholtz vortices 
(Foullon et al. 1201 It Ofman & Thompson 1201 II ) increased the 
interest in flow instability in the solar corona. On the other hand, 
the twisted magnetic field itself can cause the so-called kink in- 
stability even when there is no flow. The stability of static twisted 
magnetic flux tubes has also been studied extensively in the con- 
text of laboratory plasma (e.g., Shafranov |19571 Kruskal et al. 
11958b . Magnetic tubes are subject to the kink instability when 
the twist exceeds a critical value (Lundquist ll95U Hood & Priest 
1 1979b . Oscillations and waves and their stability in twisted mag- 
netic flux tubes without flow have been studied in the frame- 
work of the normal mode analysis in earlier works (see Dungey 
& Loughhea d[1954t Roberts [1^561 Treh an [1958l Bogdan [TM4l 
Linton et al. [19961 Bennett et al. [19991 Erdelyi & Ca rter [20061 
Erdelyi & Fe dun [20061 [20071 [201 01 Rud erman [20071 Ca rter & 
Erdelvi l2008l Ruderma n & Er delyi |2009l Karami & Barin |2009l 
and Karami & Bahari 120101 ). Most of these papers deal with 
relatively simple twisted magnetic configurations: incompress- 
ible plasma cylinders/slabs surrounded by perfectly conducting 
unmagnetised plasma or a medium with an untwisted homoge- 
neous magnetic field. Erdelyi & Fedun ( 12006b were the first to 
study the wave propagation in a twisted cylindrical magnetic 
flux tube embedded in an incompressible but also magnetically 
twisted plasma. A more complex magnetic field configuration 
has been studied by Erdelyi & Carter (2006). who investigated 
the propagation of magnetohydrodynamic (MHD) waves in a 
fully magnetically twisted flux tube consisting of a core, an an- 
nulus, and an external region. These authors considered mag- 
netic twist just in the annulus, while the internal and external 
regions retained a straight magnetic field. Two modes of sausage 
oscillations occur in this configuration, notably pure surface (i.e., 
evanescent) and hybrid (spatially oscillatory in the twisted annu- 
lus, otherwise evanescent) waves. The propagation characteris- 
tics of the kink mode in the same magnetically twisted annu- 
lus have been studied by Carter & Erdelyi (120081) . The influ- 
ence of compressibility on the dispersion characteristics of the 
sausage mode (m = 0) propagating in a magnetically twisted 
flux tube embedded in a compressible uniformly magnetised 
plasma environment in cylindrical geometry has been investi- 
gated by Erdelyi & Fedun (2007). That analysis has been gen- 
eralised by the authors (see Erdelyi & Fedun [20T0b for the kink 
mode (m = 1), as well as for the fluting modes (m > 1); m being 
the azimuthal order of a given mode. The last work by Erdelyi & 
Fedun ( 120101 ) - the most elaborate study to date - deals with the 
MHD wave propagation along twisted magnetic tubes of com- 
pressible ideal plasma surrounded by also compressible but uni- 
formly magnetised fully ionized medium. The effect of a twisted 
magnetic field on the spectra and the resonant absorption of kink 
and fluting surface modes in resistive incompressible cylindri- 
cal plasmas have been studied by Karami & Barin (2009J) and 
Karami & Bahari ( 12010b . It was found that a magnetic twist will 
increase the frequencies, damping rates, and the ratio of the os- 
cillation frequency to the damping rate of these modes. The pe- 
riod ratio P\/P2 of the fundamental and its first-overtone sur- 
face waves for kink (m = 1) and fluting (m = 2,3) modes is 
lower than two (the value for an untwisted loop) in the presence 
of a twisted magnetic field. In particular, for the kink modes, the 
magnetic twists B^/B- = 0.0065 and 0.0255 can achieve devia- 
tions from two of the same order of magnitude as in the observa- 
tions. The only work studying the wave propagation in a twisted 
magnetic tube with a mass density variation along the tube is 
that of Ruderman (2007). He investigated non-axisymmetric os- 
cillations of thin (tube radius a much smaller than the finite tube 



length L) twisted magnetic tubes in a zero-beta plasma. With an 
asymptotic analysis, Ruderman (2007) showed that the eigen- 
modes and the eigenfrequencies of the kink and fluting oscilla- 
tions are described by a classical Sturm-Liouville problem for a 
second-order ordinary differential equation. He also concluded 
that the results concerning non-axisymmetric waves in twisted 
magnetic tubes obtained by Bennett et al. (11999b for incompress- 
ible plasmas can be applied to global non-axisymmetric modes, 
i.e., kink and fluting modes, in coronal loop even though the 
coronal plasma is a low -beta plasma. This conclusion also agrees 
well with the results obtained by Erdelyi & Fedun (12007b . 

An open question is how a flow along a twisted magnetic 
flux tube will change the dispersion properties of the propagating 
modes and their stability. It turns out that the flow may decrease 
the threshold for the kink instability, as was tested experimen- 
tally in a laboratory twisted plasma column (Furno et al. 12007b . 
This observation was theoretically confirmed by Zaqarashvili et 
al. d2010b . The authors studied the influence of axial mass flows 
on the stability of an isolated twisted magnetic tube of incom- 
pressible plasma embedded in a perfectly conducting unmagne- 
tised plasma. Two main results were found. First, the axial mass 
flow reduces the threshold of the kink instability in twisted mag- 
netic tubes. Second, the twist of the magnetic field leads to the 
Kelvin-Helmholtz instability of sub-Alfvenic flows for the har- 
monics with a sufficiently large azimuthal mode number m. The 
effect is more significant for photospheric magnetic flux tubes 
than for coronal ones. In a recent paper, Diaz et al. (1201 lb also 
studied the equilibrium and stability of twisted magnetic flux 
tubes with mass flows, but for flows along the field lines. The au- 
thors focused on the stability and oscillatory modes of magnetic 
tubes with a uniform twist in a zero-beta plasma surrounded by 
a uniform cold plasma embedded in a purely longitudinal mag- 
netic filed. Regarding the equilibrium, the authors claimed that 
the only value of the flow that satisfies the equations for their 
magnetic field configuration is a super-Alfvenic one. The main 
conclusion is that the twisted tube is subject to the kink in- 
stability unless the magnetic field pitch is very high, since the 
Lundquist criterion is significantly lowered. This is caused by 
the requirement of having an Alfven-Mach number greater than 
1, so the magnetic pressure balances the magnetic field tension 
and fluid inertia. The authors suggest that this type of instability 
might be observed in some solar atmospheric structures, such as 
surges. 

Kink waves are frequently observed in the chromospheric 
spectral lines as a periodic Doppler shift and/or periodic spa- 
tial displacement of the tube axis (Kukhianidze et al. 120061 
Za qarash vili et al. [20071 Zaq arashvili & E rdelyi [20091 He et 
al. 120091 Pietarila et al. I20TT1 Tavabi et al. I20TT1 Ebadi et al. 
120121) . Most of the observations have been associated to spicules, 
which are chromospheric plasma jets. It is well known that a 
simple hydrodynamic jet is unstable against the antisymmetric 
(kink) mode (Drazin[ 2002| Zaqarashvili l201 Ik therefore the ob- 
served periodic displacements can be connected to the unstable 
kink mode of a plasma jet in spicules. However, the magnetic 
field of spicules may significantly influence the stability of the 
kink mode in the plasma jet, therefore studying the instability of 
m = 1 mode of the jet in a magnetic cylinder is of vital impor- 
tance. 

In the present work, our aim is to investigate the effect of a 
twisted magnetic field on the stability of kink (m = 1) waves 
propagating on a cylindrical incompressible radially homoge- 
neous tube with axial mass flow, assuming that the waves are 
subject, under certain conditions, to the Kelvin-Helmholtz in- 
stability. To evaluate that effect, we compare the critical flow 
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speeds at which the instability occurs with those in the untwisted 
magnetic tubes. It is necessary to point out that our simplified 
model of radially homogeneous plasma density inside the tube 
rules out the resonant absorption of MHD modes that can modify 
both the wave frequency spectrum and conditions under which 
the Kelvin-Helmholtz instability occurs. The radial density in- 
homogeneity is usually introduced as a nonuniform transition 
layer that continuously connects the internal density to the ex- 
ternal one. Owing to this transverse inhomogeneous transitional 
layer, wave modes with m + 1 are spatially damped by resonant 
absorption. In particular, the effect of longitudinal flow on the 
resonantly damped propagating kink waves has been studied by 
Terradas et al. d2010al 1201 Obi) and by Soler et al. d2010aD . We 
note that a Kelvin-Helmholtz instability of kink oscillations in 
flux tubes with both a sharp and smooth transition layer can also 
occur as a result of shear motions (Terradas et al. 120081 ). But as 
Terradas et al. (2008) claim, a magnetic twist, not included in 
their model, might decrease or even suppress the instability be- 
cause a magnetic field component along the flow stabilises the 
Kelvin-Helmholtz instability. The same conclusion for the sta- 
bilising role of a small azimuthal component of the magnetic 
field has been drawn in the paper by Soler et al. |2QlQbl). The in- 
fluence of a twisted magnetic field on the stability status of kink 
waves propagating along radially inhomogeneous solar cylindri- 
cal photospheric/chromospheric jets is beyond the scope of this 
study, and as was mentioned above, we limit ourselves to con- 
sidering radially homogeneous flowing plasmas only. 

This paper is organized as follows. In Sec. 2 we specify the 
geometry of the problem, governing equations, and the deriva- 
tion of the wave dispersion relation. Section 3 deals with the 
numerical solutions to the dispersion relation for various values 
of the input parameters and compares derived dispersion curves 
and growth rates of the unstable kink waves with those in un- 
twisted magnetic tubes. The last Sec. 4 summarises our results. 

2. Geometry, the basic MHD equations, and the 
wave dispersion relation 

We consider a magnetic flux tube with radius a and density 
Pi embedded in a uniform field environment with density p e . 
Both media (inside and outside the tube) are supposed to be in- 
compressible. The magnetic field inside the tube is helicoidal, 
Bi = (0, Bitter), Bi z (r)), while outside the magnetic field is uni- 
form and directed along the z-axis, B e = (0, 0, B e ). Both p\ and 
p e are assumed to be homogeneous. We consider the mass flow 
vq = (0, 0, vq) directed along the z-axis; thus the equilibrium 
mass flow is not field-aligned. No mass flow is present outside 
the tube, therefore the surrounding photospheric medium is con- 
sidered to be uniformly magnetised (B e z is constant), uniform 
(p e constant), and lacking mass flow at equilibrium. 

In cylindrical equilibrium the magnetic field and plasma 
pressure satisfy the equilibrium condition in the radial direction 



d 
d7 



B 



2\ 



j<t> 
pr 



(1) 



Bi${r) I rB l: (r) is a constant. Accordingly, the background mag- 
netic field is assumed to be 



B(r) 



(0,Ar,B k ) forr < a, 
(0, 0, B e ) for r > a, 



(2) 



where A, B K , and B e are constant. Then the equilibrium condi- 
tion (Q]i gives the equilibrium plasma pressure pi(r) as 



2, .2 



Pi(r) = Po ■ 



A z r 



where po is the plasma pressure at the centre of the tube. We note 
that the jump in the value of B^ir) across r — a implies a surface 
current there. 

The plasma motion is governed by the set of linearised MHD 
equations for an ideal incompressible plasma: 



d , _ _/ B {) B x 

p—V\ +p(v ■ V)l>! + V \px + 

at \ p 



1 



-(Bo ■ V)B l 



1 

— i 

A* 



(B x ■ V)B = 0, 



5* 



V x (dj x Bo) - V x (d x B^ = 0, 

V-i>i =0, 
V Bi =0. 



(3) 

(4) 

(5) 
(6) 



Here, the index '0' denotes equilibrium values of the fluid ve- 
locity and the medium magnetic field, and the index '1' their 
perturbations. Below, the sum pi + Bo ■ Bi/p in Eq. (0 will be 
replaced by pitot, which represents the total pressure perturba- 
tion. 

Assuming that all perturbations are proportional to 
g{r) exp [i {-tot + m<p + k z z)] with g(r) being just a function of 
r, and that in cylindrical coordinates the nabla operator has the 
form 

„ d „ 1 d - d „ 
V = — r + -—(f> + — z, 
or r d(p oz 

from the above set of equations one can obtain a second-order 
differential equation for the total pressure perturbation pi tot 



1 d 

+ I if + 



dr 2 r dr 



Pltot = 0, 



(7) 



as well as an expression for the radial component Vi r of the fluid 
velocity perturbation Vi in terms of pi to t and its first derivative 



11 to - k ■ Vo 



pY(aj-k- v f 



to 1 - \ dr 



Pi tot 



m \ 
Z— pitot I 



(8) 



In Eq. (0, k is the so-called wave attenuation coefficient, which 
characterises the space structure of the wave and whose squared 
magnitude is given by the expression 



Here, B;(r) = (b\, + B 1 ^ 12 = \B { \ denotes the strength of the 
equilibrium magnetic field, and p is the magnetic permeability. 
We note that in Eq. (Q]) the total (plasma plus magnetic) pres- 
sure gradient is balanced by the tension force (the right-hand 
side of Eq. ([1])) in the twisted field. We consider the special 
case of an equilibrium with uniform twist, i.e., the one for which 



where 



4A 2 a/ A 



pp [(a) - k ■ v Q ) 2 - u\ ] t 



(mA + k z B,) 



^pp 



(9) 



(10) 
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is the so-called local Alfven frequency (Bennett et al. 119991) . 
The numerical coefficients Z and Y in the expression of v\ r (see 
Eq. ([8])) are respectively 



2Aa>i 



-y^Up [(W - k ■ D ) 2 - Q) 2 A ] 



and 



1 -Z z 



As seen from the expressions for the wave attenuation co- 
efficient and the radial fluid velocity component perturbation, 
the two quantities have different values inside and outside the 
twisted flux tube owing to the different spatial structure of the 
magnetic field in both media, and to the different mass densities 
(p; and p e , respectively). It is important to note that the attenu- 
ation coefficient K e outside the tube is simply equal to k z , whilst 
that inside the tube can become a purely imaginary quantity, thus 
allowing the propagation of pseudosurface (body) waves along 
the tube. Notably the twisted magnetic field inside the tube en- 
ables the existence of these waves, which otherwise are pure sur- 
face modes (Bennett et al. 1 19991 >. 

Equation (O has two different solutions inside and outside 
the jet 

aJm{K\r) for r < a, 
a e K m (k~r) for r > a. 



Pitot(r) 



Here /„, and K,„ are the modified Bessel functions of the first and 
second kind, and a\ and a e are arbitrary constants. 

The boundary conditions have to ensure that the normal com- 
ponent of the interface perturbation 



i(o)-k- Do) 



remains continuous across the unperturbed tube boundary r — a, 
and also that the total Lagrangian pressure is conserved across 
the perturbed boundary. This leads to the conditions (Bennett et 
al. [19991) 

(11) 



(12) 



Applying boundary conditions (fTTT i and (fT2l) to our solutions 
of pitot and Vir ( an d respectively £ r ), we obtain after some al- 
gebra the dispersion relation of the normal modes propagating 
along a twisted magnetic tube with axial mass flow Do 
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Bl 1 


Pilot ?r 




Pilot _ 




pa 






pa 



[(fii-k- D ) - w Ai J jjj^ - 2nuo M 

k z aK' m (k z a) 
K m {k z a) 



Pe / , 2 \ A 2 k z aK' m (k z a) 

— [Or — CO. I H 

Pi v AeJ ppi K m (k z a) 



(13) 



This dispersion relation is a generalisation of Eq. (23) in 
the work by Bennett et al. ( 119991 ) that is applicable to a twisted 
magnetic flux tube without flow and of Eq. (13) in the work by 
Zaqarashvili et al. (2010]) that is applicable to a twisted mag- 
netic flux tube with a flow embedded in a non-magnetic en- 
vironment. As we have already seen, the wave frequency co is 
Doppler-shifted inside the jet. In the above equation a prime (/) 
denotes the derivative of a Bessel function to its dimensionless 



argument, and the local Alfven frequencies inside and outside 
the tube are correspondingly 



1 



W A i = 



(mA + k z B lz ) and a>\ e - 



k z B e 



The wave attenuation coefficient inside the tube, according to 
Eq. is given by 



Ki = k z 



4A 2 w 2 A . 

Ai 



ppi [(w - k ■ D ) 2 - a) 2 M ] / 



1/2 



Assuming no twist, when = 0, the dispersion relation 
(fT3l reduces to 



Pe ( 2 

— [(A) 



^2 \ 



(14) 



The Alfven speeds in this equation are expressed in terms of the 
background magnetic fields in both media and the corresponding 
mass densities via the well-known definitions 



t'Ai 



5, 



and 



5 e 



''Ae 



respectively. Equation (fT4l is akin to a dispersion relation in 
the same form as was previously obtained by Edwin & Roberts 
(119831 ) for untwisted magnetic tubes without flow. We also re- 
call that for the kink mode (m = 1) one defines the so-called 
kink speed (Edwin & Roberts [1983) 



Ck = 



/ 2 2 \l/2 

Pi + Pe , 



/.,2 
Ai 



1 +Pe/Pi 



1/2 



(15) 



which characterises the propagation of transverse perturbations. 

3. Numerical solutions and wave dispersion 
diagrams 

We focus our study on the propagation of the kink mode, i.e., for 
m = 1 . It is obvious that Eq. ( IT3b and Eq. (fT4l) can be solved only 
numerically. The first step is to define the input parameters that 
characterise the twisted magnetic tube and the axial flow, and to 
normalise all variables in the dispersion equations. The density 
contrast between the tube and its environment is characterised 
by the parameter rj = p e /pi- In investigating the kink instability 
in twisted magnetic tubes with axial and field-aligned flows with 
their equilibrium magnetic fields in the form B{r) = (0, Ar, Bq), 
Zaqarashvili et al. (2010) and Diaz et al. (1201 II ) chose as twist 
characteristics the so-called 'dimensionless pitch' k~p, in which 
the specific length p is defined (in their notation) as 



So 
A 



or more generally as 



~A 



That dimensionless pitch can be presented in the form k z p = 
2np/A, where 27rp is the magnetic field pitch in its common 
sense. As we showed, this parameter is associated with a fixed 
wavelength A. To study the Kelvin-Helmholtz instability in these 
tubes we prefer to specify the twist by the ratio of the two 
magnetic field components, and B z , evaluated at the inner 
boundary of the tube, r — a, i.e., via B^,/B z , where B^ = Aa. 
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This choice is symbolically equivalent to the replacement of 
the wavelength A by the tube radius a. Following Aschwanden 
d20051 l. the magnetic twist parameter B ( f ) /B z can be estimated as 
B^,/B z = tan#, where 6 is the shear angle between the untwisted 
and twisted field line. The geometric shear angle, 6, can observa- 
tionally be measured in twisted photospheric magnetic tubes and 
in twisted coronal loops. We note that our choice of the twist pa- 
rameter allows us to immediately compare the dispersion curves 
of kink waves and the growth rates when the waves become un- 
stable in twisted magnetic tubes with those in untwisted ones. 
As usual, we normalise the velocities to the Alfven speed 

UAi = Bi/inpi) 1 ! 2 , where B; = {b 2 ^ + B 2 ^j is evaluated at the 
inner boundary of the tube, r — a. Thus, we introduce the dimen- 
sionless wave phase velocity V p h = yph/^Ai and the Alfven-Mach 
number Ma = Uo/^Ai, the latter characterising the axial mass 
flow. The wavelength is normalised to the tube radius a, which 
means that the dimensionless wave number is K — k z a. It is 
worth to discuss the normalisation of the local Alfven frequency 
to a;. After multiplying both sides of Eq. (TTOb by the tube radius 
a and extracting a B\ from the brackets, we derive that 

B\ I BA B h ( B 2 / 

aa>A\ — — ■ — m vk-a — , where — = 1 r 

Vml B; ' Bi)' Bi { B\ t 

To simplify the numerical code, and largely for clarity, we can 
redefine our twist parameter as B^/Bi = e, and consequently the 
normalised local Alfven frequency is expressed in terms of the 
Alfven speed i>Ai, the dimensionless wave number K, and s: 

aa>Ai — Vai (ms + Kb z ) , 

where the dimensionless axial magnetic field component b, = 

I ,\l/2 

( 1 - e ) . The new twist parameter, e, is related to the real twist 
B^/B- via the simple relation 

B^IB Z 



s — 



^1+(B^/B z f 



(16) 



Thus, if we have an observationally measured shear angle 6 of 
some solar twisted magnetic tube, then B$/B z = tan 6 can be 
inserted into expression ( fTo*b to obtain the adopted twist s of that 
specific magnetic tube. 

Before starting the numerical calculation, we have to spec- 
ify the values of the input parameters. First we will examine, as 
in the work of Zaqarashvili et al. (120101 . the case of an isolated 
twisted magnetic tube with axial mass flow, i.e., with B e = 0, 
and hence for parameter b = B e /B\ = 0. Then, we will study the 
influence of the external magnetic field on the instability onset 
by choosing b equal to 0.1 and 0.5, respectively. Concerning the 
density contrast, rj, we consider two cases, namely 77 = 2 and 
r] = 0.1. Following Ruderman ( 120071 1. to satisfy the Shafranov- 
Kruskal stability criterion for a kink instability, we assume that 
the azimuthal component of the magnetic field is smaller than 
its axial component, i.e., we will choose our twist parameter s to 
be always less than 1 . In particular, we will study the dispersion 
diagrams of kink waves and their growth rates (when the waves 
are unstable) for three fixed values of s: (untwisted magnetic 
tube), 0.1, and 0.4. The Alfven-Mach number Ma will take val- 
ues from zero (no flow) to any reasonable number. 

3. 1 . Kink waves in twisted magnetic tubes with a density 
contrast r/ = 2 

We start by calculating the dispersion curves of kink waves as- 
suming that the twist parameter e = 0.1, and that the angular 



wave frequency w is real. As a reference, we first assume that 
the plasma in the flux tube is static, i.e., Ma = 0. The dispersion 
curves, which present the dependence of the normalised wave 
phase velocity, fph/fAi. on the normalised wave number, k z a, are 
in this case shown in Fig. Q] One can recognise two types of 
waves: a sub-Alfvenic wave labelled Ck (the green curve), and a 




Fig. 1. Dispersion curves of kink waves propagating along a 
twisted magnetic flux tube at rj = 2, B e /B\ = 0, B^/Bi = 0.1, 
and Ma = 0. The green curve is the dispersion curve of the wave 
associated with the kink speed c^. The two purple curves depict 
a band that contains nine dispersion curves of similar form, not 
plotted here. 



family of super- Alfvenic waves accommodated in a narrow band 
depicted by the two purple dispersion curves. The green disper- 
sion curve is related to the kink speed defined by Eq. (flBT) . For a 
weakly twisted magnetic tube (s = 0.1) with a density contrast 
77 = 2 the kink speed is equal to 0.5774i>Ai - our numerical code 
yields in the long wavelength limit, = 0.5738i>Ai, which agrees 
well with the value calculated from Eq. (fT3T l. If we assume that 
the Alfven speed inside the tube is typically 10 kms~', then the 
kink speed in the twisted magnetic tube is equal to 5.7 kms~', 
or approximately to 6 kilometers per second. We recall that the 
expression of Ck was defined for untwisted magnetic tubes and 
there is no guarantee that it will be valid for twisted tubes as 
well. 



Bennett et al. d 19991 1. in studying the kink wave propagation 
in a twisted magnetic tube with parameters similar to ours, had 
discovered a band that contains an infinite number of pseudosur- 
face (body) wave harmonics (Fig. 5(a) in their paper). Among 
those, generally super-Alfvenic, waves we were able to easily 
calculate only 1 1 dispersion curves from Eq. ( fT3l ) with search 
steps (in decreasing order) from 0.00025 to 0.000000025. They 
are not plotted in Fig.[T]because there is no room for all of them 
in that narrow band, but they will be shown shortly in a similar 
plot for a tube with e = 0.4. With diminishing search step, tak- 
ing it equal to 0.0000000025 or smaller, one can extract a few (3) 
additional curves, but finding them is an extremely hard numer- 
ical task - for such small search steps the data are very noisy. 
However, these high-harmonics pseudosurface kink waves are 
not too important to us since we are primarily interested in the 
evolution of the kink mode associated with the kink speed Ck 
when the plasma starts to flow. It is worth noticing that the Ck- 
labelled dispersion curve for small dimensionless wave numbers, 
till k z a = 0.1228, describes a pseudosurface (body) wave, while 



5 



I. Zhelyazkov & T. V. Zaqarashvili: Kelvin-Helmholtz instability of kink waves in twisted magnetic tubes 




< 
> 



E 




Fig. 2. (Left panel) Dispersion curves of forward propagating stable and unstable kink waves in a twisted magnetic flux tube with 
the same input parameters as in Fig.[TJfor various values of the Alfven-Mach number Ma- The non-labelled green closed curve 
corresponds to Ma = 1.325, and the black one - to Ma = 1.328. (Right panel) Growth rates of the unstable kink waves for 
M A = 1.32205, 1.347, 1.372, and 1.397. 
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Fig. 3. (Left panel) Dispersion curves of forward propagating stable and unstable kink waves in an untwisted magnetic flux tube at 
j] — 2 and B t jB\ = for various values of the Alfven-Mach number Ma- (Right panel) Growth rates of the unstable kink waves for 
M A = 1.30825, 1.325, 1.35, and 1.375. 



behind that value the wave is a pure surface mode. These types of 
waves are sometimes called hybrid waves (Bennett et al. 1 1999b . 

It is essential to emphasise that the dispersion equation ( fT3l 
also yields solutions for negative values of the wave phase ve- 
locity (i.e., for backward-propagating waves), which are a mir- 
ror image of the solutions in the positive semi-space. Including 
the flow, the whole pattern shifts en bloc upwards, changing 
of course both the Ck-labelled dispersion curve and the wave 
harmonics dispersion curves. For each Ma we have two c^- 
dispersion curves that initially, for small Alfven-Mach numbers, 
are independent, but at certain magnitudes of Ma begin to in- 
teract; for example, they form semi-closed or closed dispersion 
curves by merging. These dispersion curves are an indication 
that we are in a region in which the kink wave may become un- 
stable (see Zhelvazkov l20 1 01 ). 

To study the stability/instability status of the c^-kink wave, 
we have to assume that the wave frequency is complex, i.e., 
co — > to+iy, where y is the expected instability growth rate. Thus, 
the dispersion equation becomes complex (complex wave phase 
velocity and real wave number) and deriving the solutions to a 
transcendental complex equation is generally a hard task (see 
Acton [19901 . We numerically solve Eq. ( fT3l using the Muller 
method (Muller fl956b to find the complex roots at fixed input 
parameters 77 = 2, b — 0, and e = 0.1, and varying Alfven- 
Mach numbers Ma from zero to generally reasonable values. 
The evolution of the two kink waves with increasing flow veloc- 
ity (or equivalently Ma) is shown in Fig. [2] We note that while 
for rest plasma or small Alfven-Mach numbers the lower kink 



wave (i.e., the one that is a mirror image of the green curve 
in Fig. [TJ is a backward one, for Ma = 1 it becomes a for- 
ward wave. At the next increase of the Alfven-Mach number 
the two dispersion curves change their form and structure - ini- 
tially being hybrid modes, for Ma > 1 .23 they are pure surface 
waves. The most interesting observation is that for Ma > 1.2173 
both curves begin to merge and for Ma > 1.2933 they form 
a closed dispersion curve. The ever increasing Ma yields yet 
smaller closed dispersion curves - the smallest one depicted in 
the left panel of Fig. [2] corresponds to Ma = 1.328. All these 
dispersion curves present a stable propagation of the kink waves. 
However, for Ma > 1 .32205 we obtain a new family of wave dis- 
persion curves that correspond to an unstable wave propagation. 
We plot in Fig. [2] four curves of that kind that were calculated 
for M A = 1.32205, 1.347, 1.372, and 1.397, respectively. The 
growth rates of the unstable waves are shown in the right panel 
of the same figure. The instability that arises is of the Kelvin- 
Helmholtz type. Note the shape of red curve labelled 1 .32205: 
the growth rate has relatively high values for small normalised 
wave numbers, then passes through a minimum (its relative value 
there is ~0.003) after which it starts to grow. That specific curve 
is to some extent a marginal curve - for values of Ma lower than 
1 .32205 the kink wave is stable whilst at Ma > 1 .32205 the wave 
is definitely unstable. This means that the instability onset at the 
critical Ma will prevent the occurrence of the neutral dispersion 
curve with M A = 1.328. 

It is interesting to see how the twist of the magnetic field 
has changed the dispersion curves and growth rates of unsta- 
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ble kink waves calculated for an untwisted magnetic tube. In 
that case the dispersion curves and growth rates can be ob- 
tained through the exact numerical solutions of Eq. ( fT4b (see 
Zhelyazkov ( 120101 1). The results are shown in Fig. [3] Comparing 
Figs. [2] and [3j we immediately see that the dispersion diagrams 
and growth rates of the unstable kink waves are very similar. 
There are differences, of course; firstly the magnetic field twist 
visibly extends the closed dispersion curves in the horizontal di- 
rection, and secondly, it slightly increases the threshold Alfven- 
Mach number for the onset of the Kelvin-Helmholtz instability - 
for the untwisted magnetic tube it is equal to 1.30825, whilst for 
the twisted one it is slightly higher - the corresponding value is 
1.32205. The similarities between the dispersion diagrams and 
the growth rates are important to us - they indirectly confirm the 
correctness of our numerical solutions derived by solving the 
complex dispersion equation ( [TBI . 




o 1 1 1 1 1 1 1 1 1 
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Fig. 4. Dispersion curves of kink waves propagating along a 
twisted magnetic flux tube at r\ = 2, B e /B\ = 0, fi^/fi; = 0.4, 
and Ma = 0. The green curve is the dispersion curve of the wave 
associated with the kink speed c^. The purple curves depict all 
1 1 dispersion curves; 9 of them are divided into two sub-groups. 
The 'evolution' of those curves when increasing k z a can be seen 
in the next figure, which shows them near the left top and right 
bottom corners of the plot. 

Figure|4]shows the dispersion curves of kink waves when the 
twist parameter e = 0.4 and there is no flow. Now one can see 
all 1 1 easily derived dispersion curves of the wave harmonics. 
They have similar shapes, and as seen in Figs. |4] and [5] one can 
distinguish two sub-groups of four and five curves, respectively. 
We note that the eleventh (the highest) dispersion curve is not 
visible in the left panel of Fig. [5] because it simply begins at 
about k z a = 0.23, which is far beyond the right side of the plot. 
Another important observation is that the kink speed is markedly 
lower than expected from the Eq. (TT~5T > value - here its normalised 
magnitude is equal to 0.5 167, which corresponds to a kink speed 
of 5 kilometers per second. The kink wave itself is a hybrid mode 
- till k z a = 0.5623 it is a pseudosurface (body) wave, becoming 
a pure surface wave for the next values of the normalised wave 
number. When plasma in the tube is flowing, the kink waves are 
pure surface mode for Ma > 1.197. 

The dispersion curves and growth rates of kink waves for 
this higher value of the magnetic field twist are given in Fig. [6] 
At the left panel of this figure we see dramatic changes in the 
shape and position of the closed dispersion curves. The most 
striking changes, however, are the circumstances that the onset 




0.9 1 ' 1 ' ' 1 ' ' 1 
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Fig. 5. Zoomed-in parts of the 1 1 dispersion curves near the left 
top (top panel) and right bottom [bottom panel) corners of the 
plot in Fig [4] 

of unstable kink waves, like in Fig. [2j starts at a value of the 
Alfven-Mach number lower than that of the narrowest closed 
dispersion curve - 1.250075 vs 1.29! In untwisted flux tubes the 
kink mode appears at Alfven-Mach numbers higher than those of 
the smallest closed dispersion curves, as seen in Fig.|3](cf. also 
Fig. 4 in Zhelyazkov (2012a )). In reality, of course, the onset 
of a Kelvin-Helmholtz instability in a twisted tube will 'reject' 
the co-existence of stable kink waves like those labelled 1 .267, 
1.28, 1.287, and 1.29 in Fig. 1 

Similar calculations were performed for the two external 
magnetic fields specified by b — 0.1 and 0.5. As expected, the 
environment's magnetic field stabilises the kink wave propaga- 
tion (Bennett et al. 1999) - the instability onset occurs at slightly 
higher threshold Alvfen-Mach numbers, but the shapes of the 
dispersion and growth rate curves are similar to those of an iso- 
lated twisted flux tube as depicted in Figs. [2] [3] and [6] The in- 
fluence of the magnetic field twist, e, at a fixed external mag- 
netic field, b, and that of the environment's magnetic field at 
a fixed tube twist on the specific growth rate curves are illus- 
trated in Figs. [7] and [8] As seen in Fig. [7] at a fixed magnetic 
field configuration, the tube twist shifts to the right the mini- 
mum of the specific growth rate curve. This is most pronounced 
at b = 0.5. A twist of 0.4 notably decreases the wave growth rate 
for each b and yields the lowest threshold Alfven-Mach number 
for a given b. Thus, for a photospheric twisted magnetic tube 
with s = 0.4 flow speeds in the range 12.5-13.6 kms -1 can 
ensure the occurrence of a Kelvin-Helmholtz instability of the 
kink waves propagating along that tube. These flow speeds can 
be observed/detected in the solar photosphere. It becomes clear 
from Fig. [8] that an external magnetic field visibly increases the 
threshold Alfven-Mach number and diminishes the wave growth 
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Fig. 6. {Left panel) Dispersion curves of forward-propagating stable and unstable kink waves in a twisted magnetic flux tube with 
the same input parameters as in Fig.[4]for various values of the Alfven-Mach number Ma- (Right panel) Growth rates of the unstable 
kink waves forM A = 1.250075, 1.275, 1.3, and 1.325. 
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Fig. 7. Specific growth rates of kink waves propagating along an untwisted/twisted magnetic flux tube at 77 = 2, B e /Bi = (left 
panel), B e /B[ = 0.1 (middle panel), and B e /B[ = 0.5 (right panel) for various values of the twist parameter B^/Bi, equal to (red 
curve), 0. 1 (green curve), and 0.4 (blue curve), respectively. The curve labels denote the threshold Alfven-Mach number. 
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Fig. 8. Specific growth rates of kink waves propagating along an untwisted/twisted magnetic flux tube at 77 = 2, B^jB\ = (left 
panel), B^/Bi = 0.1 (middle panel), and B^/Bj = 0.4 (right panel) for various values of the environment's magnetic field B e /Bi, 
equal to (red curve), 0.1 (green curve), and 0.5 (blue curve), respectively. The curve labels denote the threshold Alfven-Mach 
number. 



rate - the latter effect is strongest for a magnetic tube with twist 
parameter s = 0.4. 

3.2. Kink waves in twisted magnetic tubes with a density 
contrast 77 = 0. 1 

A density contrast 77 = 0.1 is perhaps more suitable for a coronal 
and upper chromospheric magnetic flux tube than to a photo- 
spheric one, but nevertheless, studying this case is important for 
us to see how the density contrast affects the threshold/critical 
Alfven-Mach numbers for the onset of the Kelvin-Helmholtz 
instability of kink waves propagating along twisted magnetic 
tubes. The procedure for performing calculations and plotting 
dispersion diagrams is the same as in the previous subsection. 



This lower value of 77 yields markedly higher threshold Alfven- 
Mach numbers as well as higher growth rates compared to 77 = 2. 
Here, in Fig. [9] we show only one set of plots - namely that for 
an isolated flux tube (b = 0) with a twist of 0.4. At these values 
of the input parameters we obtain the smallest threshold Ma- In 
Fig. [10] like in Fig. [7] we present the specific growth rate curves 
at fixed external magnetic fields for the three different values of 
the magnetic twist. As seen from this figure, for that value of the 
density contrast (77 = 0.1) the influence of the magnetic twist, s, 
on the shapes and positions of the specific growth rate curves, ir- 
respective of the magnitude of the environment's magnetic field, 
is negligibly small. That is why a figure, similar to Fig. [8] is 
useless - in each subfigure for given s the three curves would 
practically coincide because the differences are only in the sec- 
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Fig. 9. {Left panel) Dispersion curves of forward-propagating stable and unstable kink waves in a twisted magnetic flux tube with 
T] = 0.1, B e /B\ = 0, and fi^/B; = 0.4 for various values of the Alfven-Mach number Ma- (Right panel) Growth rates of the unstable 
kink waves for M A = 3.778422,3.81,3.835, and 3.86. 
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Fig. 10. Specific growth rates of kink waves propagating along an untwisted/twisted magnetic flux tube at r\ — 0.2, B t jB\ = [left 
panel), B e /B[ = 0.1 (middle panel), and B e /B[ = 0.5 (right panel) for various values of the twist parameter B^/Bi, equal to (red 
curve), 0. 1 (green curve), and 0.4 (blue curve), respectively. The curve labels denote the threshold Alfven-Mach number. 



ond/third places behind the decimal point. As for a photospheric 
twisted tube, one obtains the lowest threshold Ma at b — and 
e = 0.4 - the corresponding value is 3.778422. If we accept 
that such a twisted magnetic tube can model a coronal loop, to 
estimate the Alfven speed inside the tube, we assume a typical 
particle number density n\ * 10 9 cirT 3 (see Ugarte-Urra et al. 
12005b and a magnetic field strength B\ = 1 1 G (see Verwichte 
et al. 120091 ); then u A i ~ 900 kms" 1 . With this high Alfven speed 
one would need a flow speed of about 3400 kms 1 for the on- 
set of an instability of the Kelvin-Helmholtz type. Since the de- 
tected flow velocities in the coronal loops are of the order of a 
few tens of kilometers per second, we can conclude that the kink 
waves in coronal loops with axial mass flows are stable against 
the Kelvin-Helmholtz instability. 

4. Conclusion 

We have studied the dispersion properties and the stability 
of the MHD kink modes running along the length of photo- 
spheric/chromospheric uniformly twisted magnetic tubes with 
axial mass flows. These were modelled as straight cylindrical 
jets of ideal incompressible plasma surrounded by a rest unmag- 
netised/magnetised also incompressible fully ionized medium. 
The wave propagation was investigated in the context of stan- 
dard ideal magnetohydrodynamics by using linearised equations 
for the perturbations of the basic quantities: pressure, fluid ve- 
locity, and wave magnetic field. The derived dispersion equation 
describes the propagation of the kink mode influenced by the 
twist of the equilibrium magnetic field B\ and the presence of 
moving plasma. The twist itself is characterised by the ratio of 
the azimuthal, B^, to the axial, B l: , background magnetic field 



components both evaluated at the inner surface of the twisted 
magnetic tube. To avoid competition from the so-called kink in- 
stability, which is largely caused by the twist of the equilibrium 
magnetic field, we chose the ratio B^/Bi- to be less than 1. For 
simplicity and clarity of the numerical code used for solving the 
wave dispersion relation (in complex variables) we used a safer 
twist parameter s defined as B^/Bi. Our s and the real twist 
B-^/B^ have generally rather close magnitudes, because they are 
related by a simple expression. The streaming plasma is char- 
acterised by its velocity Do, which is directed along the axis of 
the twisted magnetic tube. An alternative and more convenient 
way of specifying the mass flow is by defining the Alfven-Mach 
number: the ratio of the jet speed to the Alfven speed inside the 
jet, Ma = fo/fAi- The key parameters controlling the dispersion 
properties of the kink waves are the so-called density contrast, 
r\ = p e l p u defined as the ratio of the density of the environment 
to that of the twisted tube, the ratio of the two background mag- 
netic fields, b = B e /B[, the twist of the magnetic field, s, and, 
naturally, the mass flow characterised by the Alfven-Mach num- 
ber M A . 

To study the influence of the mass flow, the twist of the equi- 
librium magnetic field inside the tube, as well as the environ- 
ment's magnetic field on the dispersion properties and stability 
of the kink waves, we considered two cases with fixed density 
contrasts rj — 2 and 0.1, respectively. For a given j] and a speci- 
fied external magnetic field via the parameter b, we investigated 
the kink wave dispersion curves and the instability growth rate 
when the wave becomes unstable for two different magnetic field 
twist parameters, namely s — 0. 1 and 0.4, respectively, by gradu- 
ally changing the normalised mass flow velocity, i.e, the Alfven- 
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Mach number Ma, from zero (rest plasma) to some reasonable 
values. The flow shifts both the kink-speed curves and the dis- 
persion curves of the wave harmonics upwards. We focused our 
study on the kink-speed wave only since it is the mode that plays 
an important role in the dynamics of photospheric and chromo- 
spheric/coronal flows. 

To evaluate the effect of the adopted magnetic field twist, 
s, on the kink-speed wave dispersion curves and chiefly on the 
value of the critical Alfven-Mach number that determines that 
flow speed that ensures the onset of an instability of the Kelvin- 
Helmholtz type, we compared the wave dispersion diagrams in 
twisted tubes with those in untwisted magnetic tubes. It is well 
established that in untwisted magnetic tubes the kink waves are 
subject to the Kelvin-Helmholtz instability when Ma exceeds a 
certain critical value, depending upon the density contrast, 77, and 
on the ratio of the background magnetic fields, b, in both media 
(see Vasheghani Farahani et al. |20U9l Zhelvazkov l20T0l 120 1 2aD . 
For an isolated magnetic tube, b — 0, and, hence, the thresh- 
old value of the Alfven-Mach number depends only on the den- 
sity contrast, 77. Relatively high values of 77 yield lower threshold 
Alfven-Mach numbers while the lower values of 77 do just the 
opposite - for our chosen two 77s the corresponding critical nor- 
malised flow speeds are 1.30825 and 3.902503. The effect of the 
magnetic field twist critically depends on the value of the twist 
parameter s. This effect is strongest when e = 0.4 - then one 
observes a noticeable decrease in the threshold Alfven-Mach 
number: for a tube with 77 = 2 and b = 0, M£ = 1.250075, 
while at 77 = 0.1 it is equal to 3.778422. This substantial find- 
ing that the twisted magnetic field lowers the threshold Ma for 
the occurrence of the Kelvin-Helmholtz instability of m = 1 
mode agrees with the same statement declared in the works by 
Bennett et al. (1999) and Zaqarashvili et al. (2010). In particular, 
in twisted photospheric magnetic tubes with a density contrast 
77 — 2, magnetic field twist e = 0.4, and Alfven speed v& = 10 
kms -1 , a reliable mass flow of 12.5 kms~' can cause the onset 
of an instability of the Kelvin-Helmholtz type of m = 1 mode. 
Any none-zero environment magnetic field, B e , increases M", 
slightly for b = 0.1, and more noticeable at b — 0.5. The ob- 
servation/detection of the Kelvin-Helmholtz instability of m = 1 
mode in a coronal twisted magnetic tube (at 77 = 0.1) is gen- 
erally problematic, but a twist of 0.4 might help starting this 
instability in coronal X-ray flows, for which the parameter 77 
is slightly bigger (of the order of 0.14), and the observed flow 
speeds can be much higher, a few thousand kilometers per sec- 
ond (see Madiarska l201 II also Zhelyazkov 120 1 2bl and the dis- 
cussion therein). 

A criterion for the appearance of the Kelvin-Helmholtz insta- 
bility of kink waves in untwisted magnetic tubes is the satisfac- 
tion of an inequality suggested by Andries & Goossens (120011) . 
which in our notation reads 

M A > 1 + — — . 

For twisted magnetic tubes, and in particular for isolated ones 
(with B e = 0), this criterion provides a fairly rough prediction, 
just Ma > 1. In other words, the mass flows must be super- 
Apfvenic to trigger an instability - this seems to be a norm for 
twisted magnetic tubes with field-aligned flows (see Dfaz et al. 
1201 II ). Our numerical findings of the critical Alfven-Mach num- 
bers also agree with that rough criterion - in both cases (77 = 2 
and 0.1) the Kelvin-Helmholtz instability occurs when the axial 
mass flows are super- Alfvenic. 

On the other hand, a criterion for the appearance of the 
Kelvin-Helmholtz instability in isolated twisted magnetic tubes 



is the satisfaction of an inequality suggested by Zaqarashvili et 
al. (12010b . which in our notation is 

f7zM A > 1 + 1/77. 

For the kink mode, i.e., m = 1, this inequality gives Ma > 
1.2247 for 77 = 2 and M A > 3.3166 for 77 = 0.1. Our numeri- 
cal findings of the critical Alfven-Mach numbers agree with this 
criterion in both cases. 

The next steps in studying the Kelvin-Helmholtz instability 
of kink waves in twisted magnetic tubes are obviously (i) to in- 
vestigate the role of the transverse density inhomogeneity inside 
the tube and associated with it resonant wave absorption on the 
kink waves propagation and their stability/instability status, and 
(ii) to consider a field-aligned mass flow instead of an axial one. 
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